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ABSTRACT 


A Monte Carlo simulation of the effects of finite clouds 
on the light impulses produced by point and finite lightning 
sources within cubical, cylindrical, and spherical clouds 
shows that absorption is essentially negligible in the vis- 
ible and near infrared. The fractions of photons which 
escape various cloud surfaces are given as a function of 
position and geometry of the source. The light emission is 
high for intracloud discharges and the in-cloud portion of 
cloud-to-ground discharges because most lightning activity 
occurs at the -10 to -30°C temperature level, which, in 
turn, is usually close to the optical center of the cloud. 

The characteristic dimensions of the light escaping from a 
cloud surface are typically 60 to 70% of the cloud dimen- 
sions; and the time-broadening of an impulse by multiple 
scattering can be tens of microseconds or more. 






1. Introduction 


The optical effects of clouds have been of considerable Interest in 
the atr,ospheric sciences for many years. Most analytical and numerical 
solutions of the equation of radiative transfer have been obtained for 
situations where there are plane parallel sources of radiation Incident 
on Infinite homogeneous cloud layers. To the best of our knowledge, no 
computations have yet been published which describe the transport of 
radiation prod ced by a transient and finite light source, such as 
lightning, which Is located within finite clouds. 

Understanding the effects of clouds on lightning signals Is hnporfant 
because optical sensors have been or are currently being used to detect 
and locate discharges from ground-based stations (Kitagawa and KobayashI, 
1959; Brook and Kitagawa, I960; Krider, 1966; Clegg, 1971; Mackerras, 
1973; Kidder, 1973; Griffiths and Vonnegut, 1975), from aircraft 
(Vonnegut and Passarelli, I97S; Brook et al. , I960), and from satellites 
(Sparrow and Ney, 1968; Vorphai et at», 1970; Turman, 1977, 1978, 1979; 
Edgar, 1978, 1980; Orville and Spencer, 1979; Turman and Tettelbach, 

1980; Orville, 1981; Turman and Edgar, 1981). Clouds obviously affect 
the amplitudes of the light signals and the apparent dimensions of the 
optical source; and the multiple scattering alters the shapes of the 
light signals vs. time. 

Here we briefly describe a Monte Carlo computer simulation of the 
transport of visible and near-l nfrared light within a cloud. We give 
results of tests of the program and then calculate the cloud transmis- 
sion and absorption, the positions and angular distributions of photons 
which escape the cloud, and the time-broadening of pulses by multiple 
scattering. These calculations are done for point sources at various 
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locations within cubical, spherical, and cylindrical clouds; and the effects 
with finite channels are computed using symmetry and superposition. 


2. The Monte Carlo method 


The optical characteristics of clouds are described In works by Twomey 
St at. (1967), Rozenberg (1968), Danielson et at, (1969), Twomey (1971, 1977), 
Hansen and Travis (1974), van de Hulst (1980), and many others. Since the 
radius, a, .„f a typical cloud drop Is large compared to the wavelength, X, of 
visible light, the size parameter, 2TTa/X, is large and the scattering and 
absorption by the drops Is described by the MIe equations (van de Hulst, 

1957). Following Danielson et at, (1969) and others, we approximate the 
angular distribution of the photons which are scattered by a cloud drop by 
the Henyey-Greensteln phase function: 


p(y) = 


( 1 t 


(I) 


where p Is the cosine of the scattering angle measured In the scattering 
plane, and g, the asymmetry factor, is the intensity-weighted mean of p 
over all angles. The scattering Is assumed to be symmetric about the 
direction of the Incident photon. The probability that a photon Is absorbed 
during a scattering event is described by the single-scattering albedo, 

Sp equals unity for conservative scattering, and (I - Wp) is the probability 
that a photon Is absorbed per scattering event. 

Basically, our Monte Carlo simulation approximates the physical Inter-'^ 
actions between individual photons and the cloud drops, and is similar to 
the methods described by Plass and Kattawar (1968), Danielson et at, CI969), 
Van Blerkom C|97l), Bucher (.1973), McKee and Cox U974), Davis et at, (.1976), 
Alda (.1977), Davies CI978), and others. Photons are emitted Isotropically 
from a point source, travel a certain distance, and then interact with a 
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cloud drop at a known position within the cloud. During the interaction, 
photons are either absorbed with a probability of (I - w^) or scattered, 
with the angular distribution given by Eq. (1). This process continues 
until the photon either is absorbed or escapes from the cloud surface. 

Here we consider only interactions with the cloud particles, and we ignore 
any molecular absorption or scattering, either inside or outside the cloud. 

To derive the free paths of photons between interactions, we assume that 
the probability that a photon travels a distance x without an interaction is 

P(x) = (2) 

where A is the mean free path. If we choose a random number, r, between 0 
and I, a free path is given by 

X = ~ A£n(2’) . (3) 

To save on computing time, the scattering angles were selected at random 
from an array of size N, where each angle entry had equal probability. 

These angles were computed using Eq. (I); that is. 


N 



iL - . 9l ) dy ■■ ■ 

(I + 0 ^ - 


(4) 


or 


(I + ^^) - 




jy(i - g2)~ ^ ~ 




(5) 


2g 


so that the nth scattering angle, 6(n) is given by 

0(n) = cos“^ C(p^ + y^_|)/2] . 


( 6 ) 


For our application, N was typically 2000. When a photon escaped from the 
cloud surface, its position, direction, and total path length were stored in 
an array for subsequent output. 
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3. Cloud models 

For large size parameters, the extinction efficiency of a spherical drop 
Is approximately 2 (van de Hulst, I95’0; and the Interaction mean free path 
of photons with a uniform population of drops Is approximately: 


A 


1 

2m^N 


(7) 


where a Is the mean radius, and N Is the number density of drops. The total 
optical depth, t, of a uniform cloud layer of geometric thickness L Is 



Table I shows values of mean drop radii and number densities for a cloud with 
a liquid-water content o^ 0,42 gm/m®. The associated photon mean free paths 
and the geometrical dimensions corresponding to total optical depths of 80, 
200, and 400 are also shown. The 5 pm radius drops In Table I correspond to 
a typical continental cloud, and the larger radii correspond to maritime 
clouds. In the following, we will give the results of our calculations for 
a specific optical depth, and Eqs. (7) and (8) can be used to relate these 
data to the particular type of cloud. 


Table 1 . Geometric parameters. 


Mean- 

Radius 

Number 

Density 

N 

(cm“^) 

Mean 

Free Path 

A 


Dimension, L 
(km) 


cz 

(pm) 

i\ 

(m) 

T = 80 

T = 200 

T = 400 

5 

800 

8 

0.64 

1.6 

3.2 

10 

100 

16 

1.3 

3.2 

6.4 

15 

30 

24 

1.9 

4,8 

9.6 




Table 2 gives values of the single-scattering albedo and asymmetry factor 
for water and Ice spheres In the visible (0.45 ym) and near- infrared (0,87 ym) 
regions of the spectrum. The Wj, values for water drops were calculated from 
Mie theory using the refractive Index data of Hale and Querry (1973). The 
values for Ice particles at .87 ym were approximated using the refractive 
Index data of Irvine and Pollack (1968) at 0.95 ym, assuming the particles 
are spherical. Note that the single-scattering albedos of both water and Ice 
In the visible, where there are many emission lines In the lightning spectrum 
(Salanave, I960), are essentially unity; therefore, the scattering In this 
region Is conservative. The Sip values In the near- infrared, a region where 
there Is strong Nl radiation In the lightning spectrum and where silicon 
optical detectors have optimum spectral response (Wolfe and Zlssis, 1978), 
range from 0.9983 to 0,99997, all nearly conservative. The calculations 
which will be discussed below were performed for Sip values of 0.99990, 

0.99996, and 1.0000 with 0.99996 being regarded as typical. 




Table 2. Optical 

parameters. 



Spectral 

Radius 

a 

(ym) 

Water 


Ice 


Region 

“o 

9 


9 


5 

1 .0000 

0.85 

1.0000 

0.85 

Visible 

10 

1 .0000 

0.87 

1 .0000 

0.87 

( 0.45 ym ) 

15 

1 .0000 

0.88 

‘ 1.0000 

0.88 


5 

0.99997 

0.81 

0.99992 

0.82 

Near IR 

10 

0.99996 

0.84 

0.99990 

0.87 

( 0.87 ym ) 

15 

0.99995 

0.87 

0.99983 

0.87 


The asymmetry factors In Table 2 were taken from Hansen and Travis (1974), 
who assumed the particle size distribution discussed by Hansen (1971). The 



values range from 0,8i fo 0.80, and. In the following, we assume 0.84 Is 
typical. Since we will be concerned only with clouds thot have rather large 
optical depths, we can use the scaling or similarity principles discussed by 
van de Hu 1st and Grossman (1968), Twomey and Bohren (1980), and van de Hu 1st 
(1980, Ch. 14) to apply our results to other cases. That Is, the results 
for a given Wq, g, and t can be applied to a situation where there are other 
values, Uji', g\ and x* , provided that 

5j(l - g) = ai/x'd - g') (9) 

and 

(I - a)j,)T ® (I - Wj ') t' . (10) 

Recent work on the locations of lightning discharge channels and the 
electric charges Insld© thunderclouds suggests that most lightning activity 
Is confined to a rather narrow vertical region of the cloud near the “|0°C 
to -20°C temperature level (Teer and Few, 1974; Jacobson and Krlder, 1976; 
Taylor, 1978; Krehblel ct a7. , 1979; and Proctor, 1981). The ^l' to -20“C 
level, in turn. Is usually located just about halfway between cloud base and 
cloud top In both continental and maritime siorms (see, for example, Byers 
and Braham, 1949, and Magono, 1980). In the following, we assume that the 
light Impulses which are produced by the Individual components of an Intra- 
cloud discharge, or the In-cloud portion of a cloud-to-ground flash, are 
Isotropic point sources located on or near the horizontal mid-plane of the 
cloud. The effects of an extended channel geometry, either In the mid-plane 
or with large vertical components, can be derived from these results, using 
the symmetry of the cloud and the principle of superposition. 

4. Program verification 

The two main limitations of the Monte Carlo method are the limited number 
of photons which can be traced through an optically thick cloud with finite 
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computer time, and the fact that a Monte Carlo program cannot be formally 
proven to be mathematically correct for a finite number of photons. Here, 
we choose a sufficient number of photons to keep most of our numerical 
values accurate to at least 5/6; and we verify the program by comparing 
certain specific calculations with previous analytical and Monte Carlo 
results. In the computations to be presented below, the standard devia- 
tion of a particular numerical value Is determined using the binomial 
distribution. That Is, If the fraction of photons contributing to a 
Dart leu I ar result Is y, and If the total number of photons Is n, then 
the standard deviation Is glv<^n by 



Verification of the Monte Carlo results can be obtained by computing the 
optical properties of seml-Inf Inits cloud layers with external Illumination. 
Following Chandrasekhar (1950) and Twomey and Bohren (1980), the fractional 
absorption, a, of a layer of Isotropic scatterers can be expressed by 

a = (I - Wp)^ ^(Pp) (12) 


where 7f(pp) Is the ^?-f unction as defined by Chandrasekhar (1950) and Pp Is 
the cosine of the angle of Incidence^ Twomey and Bohren (1980) show that 
this relation also approximates anisotropic scattering (g' ^ 0) If Wp is 
replaced by a scaled single-scattering albedo derived from Eqs. (9) and (10) 


so that Eq, (II) becomes 


0 


(I - g') 5o' 

I - 


a ~ 


- t 

- ‘*'0 


I - 


f - f 

g “o 


^(bo) 




( 13 ) 


(14) 
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Figure I shows a plot of layer absorptien vs, computed using Eq. (14) and 
our Mont© Carlo results for a g' of 0.B4. The agreement Is quit© good. 

Figure 2 shows Monte Carlo calculations of the reflectance of conserva- 
tive scattering layers (5^ » I ) with the optical depths and asyrwnetry factors 
as shown. The t and g values were chosen so that the scaled optical thickness 
for conservative symmetric scattering (,g^ » 0). which Is derived from Eq. (9) 

T„ « ( I - g)’V , (15) 

8 

equals 6 In Fig. 2a and equals 12 In Fig. 2b. Note that the Monte Carlo 
results preserve the scaling relation given by Eq. (15). 

Twomey et at, (1967) used a matrix method to derive the m*=>an optical path 
length of reflected and transmitted photons vs. optical thickness for layer 
clouds. Figure 3 shows a plot of these values together with our (4onte Carlo 
results, and the agreement Is quite good. 

McKee and Cox (1974) used a Monte Carlo program to calculate the reflec- 
tance of both layer and cubical clouds with conservative scattering. Figure 4 
shows the McKee and Cox values for different optical depths, together with 
those provided by our program; and again, the agreement Is good. 

5. Results 

For this study, we have computed the optical effects of cubical, spherical, 
and cylindrical clouds with total optical depths of 80, 200, and 400. These 
computations were done for a series of point sources at the locations shown In 
Fig. 5, For the cube, the coordinates of the source points are labeled In a 
scaled Cartesian coordinate system (a;,t/,3) centered on one corner and with unity 
the length of each side. The center of the cube, for example. Is the point 
(.0.5, 0.5, 0.5). If the cube had a total optical depth of 200, the center Is 
(100, 100, 100) In optical depth units; or. If we consider the cloud with 10 pm 
radius drops In Table I, the center would be (1.6, 1.6, 1.6) In kilometers. 
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Point sources within the sphere of scaled diameter unity are taken along 
a fixed radius, starting at the center (0.5, 0.5, 0.5) and moving outward 
almost to tne outer surface (0.99, .5, .5). The cylinder has a scaled height 
and diameter of unity. Point sources within the cyllhder are labeled In 
(at, I/, a) coordinates with the center at 0.5, 0.5, 0.5). 

In order to describe the photons which escape from the cloud surface. It 
will be convenient to number the faces of the cube as shown In Fig. 5a. For 
the sphere, 6 output surfaces of equal area are defined auJ labeled as shown 
In Fig. 5b. For the cylinder, the vertical surface Is divided Into 4 equal 
sections by the « “ y end x ^ - y planes, and all surfaces are labeled as 
shown In Fig. 5c. 

The results of Monte Carlo calculations of the total cloud absorption 
and the fraction of photons which escape from the various cloud surfaces are 
summarized in the Appendix, Tables A-l, A-2, and A-3, for the point sources 
sketched In Fig. 5. The total absorption Is simply the fraction of source 
photons which Is absorbed by the cloud. The escape percentages are the 
fractions of source photons which escape through the specified surfaces; 
the sum of the total absorption and the 6 escape fractions should equal |00j{. 

The point sources sketched In Fig. 5 model localized discharge components, 
such as K-changes or channel sections, which make up the In-cloud portion of 
both cloud-1'o-g round and Intracloud lightning (Kitagawa and Kobayashl, I959; 
Kitagawa and Brook, I960; Uman, 1 969; and Brook and Ogawa, 1 977). The effects 
of clouds on extended channels can be derived from the results In Tables A-l, 
A-2, and A-3, using symmetry and superposition. Figure 6 shows sketches 
of the geometry of several extended sources that were computed using the 
cubical cloud model; the associated optical results are given In Table 3. 

The structures In Fig. 6 lie on the horizontal mid^olano; however, because 
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TABI.C 3. Computiifons of photon absorption and tscape for a cubical cloud with the finite 
source channels shown In Fig. 6. AH coaiputatlons were done for £> ■ 0.64. 


Source 

Faces I 

Escape In 
i per 

Faces 

Escape In 
X per 

Total 

Absorption 
In X 
a 

Optical Parameters 

Structure 



Cloud Face 



1 Cloud Face 


t 

A 

1.2, 3, 4 

19.0 t 1.0 

5,6 



0.99996 

80 

A 



16.2 i 0.7 


j 

■SIB, 


0.99990 

200 

A 



18.8 i 0.7 



10.7 ± 0.5 

3.6 t 0.7 

0.99996 



A 



19.S t 0.7 



11.2 ± 0.5 

0 

1.00000 



A 



17.9 ± 0.6 



8.8 ii 0.5 

11.1 t 0.5 

0.99996 

400 

At 

1.2.3.4 

17.0 ± 1.0 

5,6 

15.1 t 0.9 

1.1 i 0.3 

0.99996 

80 

f 

At 



15.5 t 0.6 



13.7 i 0.5 


0.99990 

200 

At 

*■ 


16.6 ± 0.6 



14.4 t 0.6 

■SSI 

0.99996 



At 



17.4 ± 0.6 



15.2 i 0.6 

0 




At 



15.0 i 0.6 



12.3 i 0.5 

15.6 2 0.6 

0.95996 

400 

B 

1.2.3.4 

20.3 ± 0.7 

5,6 

8.1 t 0.5 

■HM 

0.99996 

200 

Bt 

1 


17.7 t 1.0 



12.5 i 0.9 

wBVi iL'^^B 

maam 

0.99996 

• 



C 

1. 

.2 

26.4 i 1.1 

7,4,5,6 

11.6 ± 0.6 

0.8 i 0.2 

0.99996 

80 

C 



26.7 t 0.7 



9.5 ± 0.5 

8.5 i 0.5 

0.99990 

1 

200 

C 



26.6 ± 0.7 



10.7 t 0.5 

3.6 t 0.3 

0,99996 



C 



27.7 t 0.8 



11.2 ± 0.5 

0 

1.00000 



C 



27.1 t 0.8 



8.8 ± 0.5 


0.99996 

400 . . 

• Ct 

1, 

,2 • 

19.0 t 1.Q 

3,4,5,6 

15.1 t 0.9 


0.99996 

80 

ct 



17.4 i 0.6 



13,7 i 0.5 


0.99990 

200 

ct 



16.6 t 0.7 



14.4 ± 0.6 


0.99996 



ct - 



19.6 t 0.7 


1 

15.2 1 0.6 

0 

1.00000 



ct 



17.7 t 0.6 



12.3 t 0.5 

15.6 t 0.6 

0.99996 

400 

D 

1, 2,3,4 

20.3 i 0.7 

5,6 . 

6.1 t 0.5 

2.7 i 0.3 

0.99996 

200 

Dt 


1 

i 

17.7 t 1.0 


1 

12.5 i 0.9 

4.2 ± 0.5 

0.99996 


1 































of tho symmetry, the results for vertical channels can be determined simply 
by relabeling the output faces. Atnong the things ts note In Table 3 are 
the comparatively small values of total absorption, even for clouds with an 
optical thickness of 400, This Is because the single-scattering albedos In 
Table 2 are all close to unity; therefore, even under extreme circumstances, 
the visible and near-infrared light escaping from the cloud will always be a 
sizable fraction of the total light emitted by the discharge. Figure 7 shows 
the total absorption vs. point source position for a nonconservative cubical 
cloud with optical depths of 200 and 4Q0. Total absorption is the greatest 
for point sources near the center of the cloud (also see Tables A-l, A-2, and 
A-3) because of the exponential behavior of the scattering process. 

Since absorption Is essentially negligible, the primary optical effects 
of clouds are a redistribution of the numbers and directions of photons which 
escape the various surfaces and the tlme=broad©nlng of the pulses due to 
multiple scattering. The fractions of photons which escape face 2 In the 
cubical cloud are plotted In Fig. 8 as a function of source position, single- 
scattering albedo, and optical depth. The nearly exponential Increase of the 
escape fraction as the source approaches the surface Is not affected by changes 
In the single-scattering albedo. An Increase In the absorption or optical 
depth merely tends to lower the entire curve. It should be noted that, w'th 
the source essentially at the surface of face 2 (0,99, 0.5, 0.5) and a t of 
200, almost 93^ of the photons escape from that face and only about 0.3^ 
escape from the opposite surface (Face I), The escape from face I for other 
source positions can be obtained simply by relabeling the a-axis In Fig. 8 to 
be from 1.0 to 0 rafher than from 0 to 1.0. The data for escape as a function 
of source position along the x - y diagonal (z = 0.5) ore also plotted In 
Fig. 8. 


KKXT--g.-j-.c- 
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For the sources plotted in Fig. 8 along the {x, 0.5, 0.5) axis, faces 
I and 2 are normal to the source axis. Faces 2, 4, 5, and 6 are parallel 
to this axis. Figure 9 shows the sum of the fractions of photons which 
escape through the parallel faces when the source moves along the 
(x, 0.5, 0.5) axis. The curves are symmetric about x = 0.5; and, again, 
the effect of increasing absorption or optical depth is io simply lower 
the entire curve. 

One of the greatest difficulties with the Monte Carlo method is in 
obtaining reliable statistics for the angular distributions of the escaping 
photons because of the large number of "boxes” into which the photon angles 
are divided. Fortunately, the nature of the statistical error is such that 
it iS smallest where the frequency of emission Is largest [see Eq. CM)]. 
Figure 10 shows the zenith angle distribution of photons which escape 
through parallel faces 5 and 6 (the bottom and top, respectively). Half of 
the escaping photons emerge within about 45 degrees of the vertical axis 
£|cos(^).| ^ .7], and the zenith-angle distributions are essentially inde- 
pendent of the source position and total optical depth. The distribution 
of azimuth angles of the photons which escape from faces 5 and 6 with the 
source positions of Fig. 10 are essentially isotropic. 

Figures |i and 12 show the zenith and azimuth angle distributions for 
photons which escape from faces I, 2, 3, and 4. The solid curves in 
Figures 10, II and 12 have actually been drawn to provide the best fit to 
the means of the angular distributions derived from all point-source loca- 
tions that are at least 0.25 scaled-distance units away from the surface. 

The enhancement of photon optical paths within the cloud by multiple 
scattering will cause a time broadening of an impulsive light signal. To 
determine the magnitude of this effect. Figures 13 and 14 show distributions 
of the optical paths of all photons which escape from the cloud surface in 
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units of the mean free path (Eq. 7). If the point sour':e Is a delta function 
with time, these curves show the time broadening when the optical path Is con- 
verted to geometric units and divided by the speed of light. That Is, each 
unit of optical path corresponds to a time Interval given by 

= A/o (16) 

where A Is given by Eq. (7) and q Is the speed of light. Since A for the 
10 ym radius drops In Table I Is 16 m, In this case would be about 

53 nanoseconds per optical path unit. 

Figure 15 shows the mean and the most probable optical paths that were 
computed for photons which escape the cubical cloud as a function of the 
poln -source location along the (a, 0.5, 0.5) axis. It should be noted In 
Figs. 13, 14, and 15 that the time-broadening Is largest when the sources 
are near the middle of tho cloud and when the optical depths are large. 

For a point source at the center (0.5, 0.5, 0,5), the mean optical paths 
of photons reaching the surface would be 970 and 2790 for cloud depths of 
200 and 400, respectively. For the 10 ym radius drops in Table I, these 
mean optical paths would correspond to moan time-delays of 51 and 148 
microseconds, respectively. It should be noted that the curves in Fig. 15 
will be symmetric about the x = 0.5 position. 

The number of photons which escape a unit area of cloud surface Is a 
function of position on the surface. In order to describe this function for 
the different faces of our cubical cloud, we will consider the dimension, 5, 
of a hypothetical square centered on each face and examine the fraction of 
photons emitted from within Its boundaries. For the spherical cloud, we 
Imagine a surface hounded by a circle of diameter, 5, on the hemisphere 
containing the source and determine the fraction of the photons which escape 




14 


this surface. A point source at the center of the sphere produces a uniform 
Illumination of the entire surface. Figure IC shows the fractions of photons 
which escape from surface 2 of the cube and the source hemisphere as a func- 
tion of the characteristic distance, S, for various point sources within 
cubical and spherical clouds. All curves lie between the cases where there 
Is uniform Illumination (curve 9) and the image of a point source on the 
surface (curves I and 2). 

If we define a characteristic dimension of the emission from a surface, 

S^, to be that value of S which contains 70^ of the photons emitted from the 
face of a cube, or, in the case of a spherical cloud, the S which contains 
10% of the photons which escape from the hemisphere containing the source, 
then we can plot as a function of the source position as shown In Fig. 17. 
Note that the characteristic dimensions of the surface Illumination are all 
about 0.6 to 0.7 scale units when the sources are deeper than about 0.2 from 
the surface. It should also be noted that the choracterlstic dimens'-Ion 
seems to be inversely proportional to the fraction of photons which escape 
through the surface (see Table A-l), 

If the cloud Is observed from a large distance, such as from a satelVite, 
the number of photons which escape from the cloud as a whole determine the 
brightness of the source. In this case, each face of the cube can be approxi- 
mated by a point source emitting the fraction of photons listed In Table A-l, 
with the angular distributions plotted In Figs, 10, 1 1 , -and 12. The fraction 

of photons, P„, passing through a certain solid angle, Ao). Is obtained by 
I 

summing over the visible faces, 

where is the fraction of photons emitted from face n through Aw^. 


u, 


n 


n 


Ao) . 


( 17 ) 
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F^(Aw.) can be expressed as F„*0 where is the fraction of 

source photons emitted from face n, is the fraction of photons emitter 

from face n with p. = cos6. between u.-^Au and u. t %Au and d) (A.) Is the 

1. t/ K t 

fraction of photons emitted from face n with 4> between and 4>. + 

Here we let Ay = .1 and A<J> = .Itt, so there are 400 solid-angle elements each 
of .OItt steradlans. 

Figures IB and 19 give the fractions of source photons which would be 
detected In the y = 0,5 plane as a function of the zenith angle, 2, measured 
from the positive s-axls and moving initially toward the positive a;-axis. 

Note that there Is a uniform brightness with angle when the source is in the 
center of either the cube or the sphere. The light intensity changes very 
rapidly with angle when Z is near 0 or tt with both cubical and sphericEil 
geometries. 

Figure 20 shews tiie angular distributions of source brightness with 
several extended channel sources, as sketched In Fig. 6. Note In Fig. 20 
how the effects of a large source geometry tend to eliminate the rapid 
intensity variations with angle. If the zenith angle, 2, is defined so that 
It moves toward the positive j/-axis, rather than the a;-axis, the brightness 
distributions for the sources plotted in Figs. 18, 19, and 20 are uniform 
with angle. 

6. Discussion 

In order to perform the above calculations, we have been forced to make 
several simplifying assumptions about the geometry and composition of thunder- 
storm clouds and the nature of lightning sources. Spheres, cubes, and 
cylinders, of course, do not exhibit the ebullient irregularities of natural 
clouds, A homogeneous composition virtually ignores any possible differences 
In the optical characteristics of water and Ice particles and possible varia- 
tions In drop concentrations and mean radii within the cloud. The 
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Henyey-Greenstein phase function Is a good approximation for most drop-size 
distributions, but it is not exact. In spite of these diff icuities, how- 
ever, we expect the basic aspects and trends of most of our results to be 
val Id. 

By modeling lightning as a series of point sources within the cloud, we 
obviously Ignore any contributions which would arise from channels ouTSide 
tho cloud surface. Light from a cloud-to-ground channel, for example, could 
reach an optical detector directly or by reflections from (or transmissions 
through) cloud surfaces. In Fig. 2 and Table 3, we see that the main effect 
of even a thin cloud. layer on the light produced by a source outside the 
surface is reflection. Table 3 also shows that the transmission through 
even a small cloud (t = 80) is limited to just a few percent of the source 
or less. 

We will now summarize our results by considering the problem of detecting 
lightning from an earth-orbiting satellite with an optical sensor. The high 
reflectivity of even thin cloud layers will prevent almost all the light pro- 
duced below or behind a cloud from reaching the detector; therefore, we can 
expect a very poor detection efficiency on most cloud-to-ground channels. 

This situation might be alleviated somewhat by reflections of light from the 
ground or other clouds, but, in these cases, the location accuracy will be 
poor. 

As we have said oefore, most lightning activity occurs within the cloud 
at altitudes which are close to the optical center. Most cloud-to-ground 
flashes are initiated inside the cloud at these same altitudes and have 
extensive in-cloud components. Since the effects of absorption tend to he 
negligible, any channels or portions of channels which are symmetric about 
the center of the cloud will produce time- integrated signals which are 
essentially the same as if there were no cloud present at all. if the 
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channels are above the geometrical midplane, the effects of the cloud will 
be to enhance the signal detected by the satellite, and sources below the 
midplane will be reduced In Intensity quite rapidly with distance. Abso- 
lute measurements of the source Intenstly will not be possible unless the 
locations of the channels Inside the cloud are known with good accuracy. 

For deep sources, the characteristic dimensions of the light escaping the 
cloud will be roughly 60 to 10% of the cloud dir.ensions; and If the channels 
are extensive, the Intensity will be substantially Independent of direction. 

The rlsetlmes of lightning light signals are probably dominated by the 
geometrical growth of the channels which produce the emissions (Krlder, 1966), 
Short channels can produce fast rlsetlmes (Krlder, 1974), but our calculations 
suggest that these rlsetlmes cannot be measured outside the cloud because the 
1 1/ne-broadenlng due to multiple scattering can easily be several tens of micro- 
seconds. 
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APPENDIX 


TAt .K A-1 . ConpuUtlons of photon ibtorptlisn and tscapt for cubital clouds 


16.4 t 0.9 
B.7 jtO.7 

6.2 t 

3.4 1 

1.3 i 0.3 


i|Cj 


farcant Escape Through Cloud lUca 


2 


16.4 tO.9 

30.7 i1.2 

39.1 i 1.3 

57.5 ±1.3 
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Tamlk A-2. CompuUtilons Of photon absorption and tscapa for a spherical cloud. 


Point 

Source 

Location 

Percent Escape Through Cloud Face 

ToUl 

Absorption 
in i 
a 

Optical Parameters 

1 

2 

3 

4 

5 

6 

h) 

• 

9 

r 

•SOtoSOfoSO 

16.1 4 0.5 

16.1 4 0.5 

16.1 4 0.5 

16.1 4 0.5 

16.1 4 0.5 

16.1 4 0 5 

3.4 4 0.3 

0.9! 

)996 

O.C 

14 

2( 

K) 

.6S,.SO,.SO 

7.2 4 0.4 

31.9 4 0.7 

14.5 4 0.5 

14.5 4 0.5 

14.5 4 0.5 

14.5 4 0.5 

3.1 4 0.2 







.7S,.S0,.50 

4.2 4 0.3 

47.8 4 0.7 

11.3 4 0.4 

11,3 4 0.4 

11.3 4 0.4 

11.3 4 0.4 

2.6 t 0.2 







.8S,.5O,.S0 

2.0 4 0.2 

69.5 4 0.7 

6.6 4 0.4 

6.6 4 0.4 
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6*6 t 0a4 

2.1 4 0.2 







.99..S0..50 

0.8 4 0.2 

95.3 4 0.7 

0.9 4 0.2 
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0.9 4 0.2 
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TAite A-3. Computations of photon absorption and escape for a cylindrical cloud. 
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FIGURE CAPTIONS 
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Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Fig. 6 
Fig. 7 
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The solid curve shows +he fractional absorption, a, of a semi- 
infinite cloud layer with vertical incidence (Tworney and Bohren, 
1980), The dots and bars a»’e the means and standard deviations 
computed using our Monte Carlo simulation. 

a) Distribution of reflected photons vs. zenith angle, Z, from 
conservative cloud layers with normal incidence. The fraction Is 
for a cos(2) Interval of 0.1. The scaled optical depth In all 
three examples equals 6. 

b) Same as a) except the scaled optical depth is 12, 

The mean optical path length, f, of reflected and transmitted 
photons vs,» optical depth from Tworney et aZ. (1967) [solid curves] 
and the results of Monte Carlo simulations [solid symbols]. 

Reflectance of vertically incident photons as determined by the 
Monte Carlo program of McKee and Cox (1974) [solid and dashed 
curves] and our Monte Carlo results [solid dots]. 

Sketches of (a) the cubical, (b) the spherical, and (c) the 
cylindrical cloud geometries used in this study. The numbers 
label the cloud faces, and the solid dots show the point-source 
locations computed for that geometry, 

(a-h) Top views of the channel geometries derived from point 
sources in the horizontal plane (s = 0.5). 

Total absorption vs. position for point sources located along 
the line (a?, 0,5, 0.5) in a cubical cloud. Note that the 
absorption is symmetric about x = 0,5. 
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Ftg. 8. The fraction of photons escaping face 2 of the cube vs. point 
source position along the ix, 0.5, 0.5) axis. 

Fig. 9. The sum of the fractions escaping through faces 3, 4, 5, and 6 
of the cube vs. source position along the line (x, 0.5, 0.5). 

Fig. 10. Distributions of emerging zenith angles, Z, for photons which 
escape through faces 5 and 6 of the cube. The fractions are 
for cos(K) Intervals of 0.1. 

Fig. II. Same as Fig. 10 except for faces I, 2, 3, and 4, 

Fig. 12. Distributions of tfie azimuth angles of photons which escape 

through faces I, 2, 3 and 4. The fractions are for azimuth I 

intervals of O.Itt radians. 

I 

Fig. 13. Distributions of the optical paths traveled by photons emitted 

from various point sources within a cubical cloud of total . i 

optical depth 200. 

Fig. 14. Same as Fig, 13 except for a total optical depth 400. 

Fig. 15, The most probable optical path, T^, and the mean optical path t 

for photons emitted from point sources along the line (x, 0.5, 0.5) 

In cubical clouds of the specified optical depth. 

Fig. 16. The fraction of photons which escape through a region of size S, 

on face 2 of a cubical cloud; and the hemispheric fraction of | 

photons which escape through a circular region of minimum diameter, I 

i 

S, on the source hemisphere. All sources are located along the | 

(x, 0.5, 0,5) axis; Cu = cube and Sp = sphere, 1 
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Fig. 17. The characteristic dimension, S^, of the region which contains 
70% of the photons emitted from cubical and spherical clouds as 
a function of point-source position along the line ix, 0.5, 0.5). 

Fig, 18. Angular dependence of the total cloud Intensity for a cubical 

cloud In the y = 0.5 plane. The fractions are for a solid angle 
Interval of O.Olir steradlans at cos(2) Intervals of 0.1. 

Fig. 19. Same as Fig. 18, except for a spherical cloud. 

Fig. 20. Same as Fig. 18, except for the extended channel sources sketched 
In Fig. 6. 
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